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Earth's proton aurora occurs over a broad MLT region and is produced by the precipitation of low-energy 
(2-10 keV) plasmasheet protons. Proton precipitation can alter chemical compositions of the atmosphere, 
linking solar activity with global climate variability. Previous studies proposed that electromagnetic ion 
cyclotron waves can resonate with protons, producing proton scattering precipitation. A long-outstanding 
question still remains whether there is another mechanism responsible for the proton aurora. Here, by 
performing satellite data analysis and diffusion equation calculations, we show that fast magnetosonic waves 
can produce trapped proton scattering that yields proton aurora. This provides a new insight into the 
mechanism of proton aurora. Furthermore, a ray-tracing study demonstrates that magnetosonic wave 
propagates over a broad MLT region, consistent with the global distribution of proton aurora. 

Earth's proton aurora is formed when charged protons precipitate into the atmosphere loss cone, within a few 
degrees at the equator, and subsequently collide with the neutral atmosphere at low altitudes \ Proton 
aurora can provide important information for understanding magnetosphere-ionosphere interaction by 
supplementing the direct imaging of the magnetosphere^. Since charged protons are easily trapped inside the 
Earth's magnetic field due to a minimum magnetic field existing at the equator, such proton precipitation requires 
a scattering mechanism to break the first adiabatic invariant. In the basically collisionless and tenuous magneto- 
sphere, interactions between protons and plasma waves can induce such non- adiabatic scattering. One important 
plasma wave, electromagnetic ion cyclotron (EMIC) wave can interact with protons and efficiently scatter 
protons into the atmosphere^"^. However, strong EMIC waves are present primarily along the plasma plume 
in the duskside^"^ and on the dayside in the outer magnetosphere L > 6^°"^^. This appears to be difficult for EMIC 
wave alone to explain the broad MLT distribution of proton auroral emission from the morning sector to the dusk 
sector associated with the precipitating protons at lower L-shells. Hence, the fundamental and long- outstanding 
problem as to what mechanism is primarily responsible for proton auroral emission still remains unresolved. 
Another important wave, fast magnetosonic (MS) wave, also named equatorial noise^^ can resonate with protons, 
potentially leading to rapid pitch angle scattering of protons. Furthermore, MS waves can propagate eastward 
(later MLT) or westward (earlier MLT) over a broad region of MLT^^"^^. However, it has not been possible so far to 
determine whether MS wave is indeed another generating mechanism of the proton aurora, because simultaneous 
observations concerning MS wave activity, proton pitch distribution and proton auroral emissions are challen- 
ging due to extremely difficult observational conditions. Fortunately, such simultaneous observations were 
serendipitously found in the unique events on September 16, 2003, to identif)^ such mechanism. 



Results 

On September 16, 2003 an interplanetary coronal mass ejection (CME) interacted with the Earth's magneto- 
sphere, leading to a small geomagnetic storm (D^^ ~ — 50 nT). The Cluster spacecrafts, which stayed around the 
noon near the Earth's equator for about ten minutes then, recorded strong MS wave activities for a duration of ~ 
10 minutes in the frequency range ~ 10-150 Hz and around the noon local time sector from 1 1.47 to 1 1.41 MLT 
(Fig. la). MS wave is confined within a few degrees of the geomagnetic equator and exhibits a notable harmonic 
structure, viz., spaced at 4/^^ and 12/^^ (Fig. la). Moreover, MS wave has a normal angle ^ ~ 90° (Fig. lb) and an 
ellipticity ~ 0 (Fig. Ic), implying that MS wave propagates very obliquely with respect to the ambient magnetic 
field direction 
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Figure 1 | Satellite data on 16 September 2003 storm. Data collected by the Cluster STAFF instrument during —10 minutes duration for MS wave power 
(a), wave normal angle (b), the angle between the Earth's magnetic field and the normal to the plane of the wave; and wave ellipticity (c), the degree of 
elliptical polarization. MS waves are right hand polarized electromagnetic waves which occur as a series of narrow tones, spaced at multiples of the 
proton gyrofrequency/cp up to the lower hybrid resonance frequency /^^h- a, MS waves maximizes basically at frequencies from ~ 10-150 Hz, spaced at 4/^^ 
ifcp ~ 5 Hz) and 12/^^. The white dotted line denotes the lower hybrid resonance frequency y^^- b-c. The observed waves have a high normal angle 6 ~ 90° 
and a high degree of elliptical polarization (ellipticity ~ 0), indicating that the k vector is approximately perpendicular to the ambient magnetic field 
direction, d. Pitch angle distribution of protons for different indicated energies (2-10 keV) measured by CIS instrument. Proton fluxes peak at a pitch 
angle of 90° and drop dramatically at small pitch angles, e-j. Auroral snapshots for northern hemisphere and southern hemisphere as a function of 
magnetic local time (MLT) and magnetic latitude (MLat) obtained by FUV-SI12 onboard IMAGE spacecraft when IMAGE spacecraft travels at the 
location ~ 15 Mlat and altitude=8 Re. Proton aurora bands are present from the morning sector to dusk sector 09:00-18:00 MLT, with the strongest 
emission basically in 14:00-18:00 MLT. The auroral emission is asymmetric with a more intensity and a broader latitude in the northern hemisphere than 
those in the southern hemisphere. 
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Figure 2 | Ray-tracing of MS waves. Ray paths representing MS waves are initiated at the equator, at MLT= 12:00 (a) and 15:00 (b), with the ray-tracing 
parameters as shown in Table 1. The corresponding color scale indicates the initial wave frequency of each ray in Hz. The black circles indicate the 
initial positions L = 5.0, 5.5 and 6.0, where L is the distance in Earth radii {IRe = 6, 370 km) from the centre of the Earth to the equatorial crossing of a 
given magnetic field line. MS waves can propagate eastward (later MLT) or westward (earlier MLT) basically from the morning sector to the dusk sector 
0900-1800 MLT, corresponding to the occurrence of proton auroral emission. 



Figure Id shows a pancake distribution of proton flux for different 
indicated energies (2-10 keV) measured by CIS instrument. Proton 
flux maximizes at a pitch angle of 90° and drops remarkably at small 
pitch angles at each energy, particularly for energy above 3.53 keV. 
Since a pancake distribution is produced when protons with smaller 
pitch angles have been scattered into the loss cone, any endeavor 
to identify whether MS waves are really responsible for proton 
auroral precipitation must also explain the formation of pancake 
distributions. 

Figure le-j show auroral snapshots observed by FUV-SI12 instru- 
ment onboard IMAGE spacecrafts^ for the northern hemisphere and 
southern hemisphere produced by the precipitation of 2-10 keV 
protons originating from the plasma sheet. The SI 12 instrument 
detects the Doppler-shifted Lyman- a photons corresponding to pre- 
cipitating charge-exchanged protons with energies of a few keV^^'^". 
Proton aurora covers a broad MLT region 09:00-18:00 MLT, with 
stronger emissions roughly in 14:00-18:00 MLT. Such a broad MLT 
distribution of proton aurora requires MS waves which can produce 
proton scattering to occur over a similar broad MLT region. 

Unfortunately, there is no direct wave data over such a broad MLT 
region during this event. Here, we adopt the previously developed 
programs^S'^^ to trace MS waves with different frequencies based on 
the wave data (Fig. la-b). MS waves are launched at 12:00 and 15:00 
MLT for different locations at the geomagnetic equator. The mod- 
eled results (Fig. 2) confirm that MS waves can indeed propagate over 
the similar MLT region from the morning sector to the dusk sector 
09:00-18:00 MLT. 

Resonant interactions between MS waves and protons take place 
when the wave frequency equals a multiple of the proton gyrofre- 
quency in the proton reference frame. These cyclotron MS-proton 
interactions produce efficient exchange of energy and momentum 
between waves and protons, leading to stochastic proton scattering. 
Such stochastic scattering can be evaluated in terms of pitch angle 
and momentum diffusion coefficients induced by cyclotron MS-pro- 
ton interactions. Then, the dynamic evolution of the proton distri- 
bution function can be calculated by solving the Fokker-Planck 
diffusion equation associated with pitch-angle and momentum 



diffusion coefficients^^. Here, we use this method to model evolutions 
of the proton distribution following their injection into the 
magnetosphere. 

Calculation of the diffusion coefficients requires a detailed know- 
ledge of the amplitudes and spectral properties of the waves. A stand- 
ard way is to assume that the wave spectral density obey a Gaussian 
distribution in wave frequency and wave normal angle^^ (see the 
details in Methods). To allow the data fitting more reasonable, we 
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Figure 3 | Gaussian Fitting curves. The modeled Gaussian fit (solid) to the 
observed two-band wave spectra (dot) over a three-minute period 
17:58:25-18:01:25 is shown, together with the fitted wave amplitudes B^, 
the peak wave frequencies the bandwidths Sfy the lower and upper 
bands fi and /2. 
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Figure 4 | Diffusion rates Proton bounce-averaged pitch angle a), momentum {{Dp^, b) and cross ((i)«p), c) diffusion rates for resonant MS 

wave interactions with protons. Pitch angle (Do^a) and cross (D^p) are respectively higher than momentum {Dp^. Combined scattering by all three 
diffusion rates leads to rapid pitch-angle scattering between 2 and 10 keV at lower pitch angles and the resultant proton auroral precipitation into the 
atmosphere on a timescale comparable to tens of seconds. Evolution of the proton flux Starting with an initial condition representative of protons ( d) , we 
show the proton flux due to scattering by MS waves from a numerical solution to the two-dimensional Fokker-Planck diffusion equation. MS wave causes 
a rapid loss of low-energy (2-10 keV) protons within tens of seconds (e-f), leading to the strongest proton auroral precipitation from the morning sector 
to the dusk sector 0900-1800 MLT. 



average the observed wave magnetic field intensity over the indicated 
time period in this event and then present a least squares Gaussian fit 
to the observed two-band spectral intensity (Fig. la), together with 
the corresponding fitting parameters as shown in Figure 3. We then 
calculate bounce-averaged diffusion coefficients for MS waves as 
shown in Figure 4a-c. All the diffusion coefficients cover a broad 
region of pitch angle and energy, with high values at lower pitch 
angles. In particular, the pitch angle diffusion rate above ~ 7 keV 
can approach 0.0 1/s around the loss-cone, allowing efficient pitch 
scattering of protons into the loss-cone in a time scale of tens of 
seconds. Moreover, pitch angle and cross diffusion coefficients are 
higher than momentum diffusion coefficients particularly at lower 
pitch angles, implying that pitch angle diffusion dominates over the 
energy diffusion while cross diffusion coefficients should also con- 
tribute to proton-MS interaction. 

Using those diffusion rates in Figure 4a-c, we solve a two-dimen- 
sional Fokker-Planck diffusion equation by a recently introduced 
hybrid difference method^^ and calculate phase space density/evolu- 
tion of protons due to MS waves. We then simulate the evolution of 
differential flux j by the subsequent conversion; = p^and extend the 
results to the range of 90°- 180° due to the mirror symmetry. The 
results (Fig. 4d-f) confirm that MS waves can produce a rapid pitch 
angle scattering of protons within tens of seconds and substantially 
modify the entire population of injected protons, leading to the 
formation of pancake distributions at energies 2-10 keV. It should 
be mentioned that anomalous depressions (or dips) occur around 



90° because of the absence of resonant scattering at (or around) 90°. 
Moreover, the diffusion rates around 90° are very small in the energy 
range 1-3 keV and gradually increase above 3 keV. Consequently, 
those dips around the pitch angle 90° are more noticeable when 
energy increases (Fig. 4a-c). Although our simulation results are 
sensitively dependent on adopted model parameters and display 
anomalous depressions around 90°, the dominant features of 
Figure 4f are very similar to the Cluster data shown in Figure Id. 
Considering that MS waves can propagate over a broad MLT region 
in this event (Fig. 2), we therefore conclude that MS waves are indeed 
responsible for the precipitation of energetic protons into the atmo- 
sphere, producing the resultant proton aurora (Fig. le-j). 

Discussion 

The simultaneous observation and corresponding modeling pre- 
sented in this study firstly link MS waves to the origin of proton 
aurora. Our results definitively demonstrate that, as presented in 
Fig. 5, MS waves can produce rapid proton precipitation responsible 
for proton aurora, naturally accounting for the remnant pancake 
proton distribution left behind in space and the broad MLT distri- 
bution of proton aurora. Although our simulations were performed 
around the noon 1 1:30 MLT, this conclusion should be valid over the 
entire region of excited MS waves since the basic properties of MS- 
proton interaction should not change no matter whether the ampli- 
tude and morphology of MS waves are different under different 
locations. 
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Figure 5 | Protons scattered by MS waves. Schematic diagram showing the solar wind protons, MS waves and the proton aurora. The image of the Earth 
was created by HZG and XFL. (1) The energetic protons originating from the plasmasheet are trapped in the Earth's magnetic field and drift westward 
around the Earth. (2) MS waves propagate eastward or westward over a broad region from the morning side to the dusk side. (3) Wave-particle 
interaction occurs in the same region. (4) It acts as an unexpected mechanism to scatter the trapped protons into the atmosphere, yielding the proton 
aurora. 



As shown in Figure 1, the proton auroral emission maximizes in 
the afternoon sector, which appears to be inconsistent with the recent 
statistical survey of MS waves observed on THEMIS spacecrafts^ that 
such waves are strongest in the dawn to pre-noon sector. However, 
the proton aurora emission intensity primarily depends on two fac- 
tors: the number of the precipitating protons and the MS wave 
intensity. Since the precipitating protons originating from the 
plasma sheet drift westwards around the Earth (see Fig. 5), the pro- 
tons encounter resonance with the afternoon sector MS waves at first, 
allowing part of protons precipitating into the atmosphere. Then the 
rest of protons continue resonating with MS waves on other sectors 
in their drifting. Moreover, the potentially existing plume EMIC 
waves along the duskside (though not observed directly here) may 
also contribute to the proton scattering. This probably explains why 
the proton aurora peak in the afternoon sector, instead of in the dawn 
to pre-noon sector. 

It should be point out that, we use the quasi-linear theory of wave- 
particle interaction. Analogous to the Particle-in-cell treatment, the 
quasi-linear theory has been frequently adopted by the space plasma 
physics and magnetosphere research community to treat wave-par- 
ticle interaction. The quasi-linear theory is valid by assuming that 
each individual particle randomly walks in velocity space, resonates 
with a succession of uncorrelated waves, and is scattered in a small 
amount in pitch angle and energy each time. Those conditions are 
basically satisfied in the Earth's radiation belts for naturally gener- 
ated MS waves, where the bandwidth is generally above the proton 
cyclotron frequency up to the lower hybrid frequency. Moreover, this 
study is intended to propose a new mechanism of proton aurora by 
MS-driven scattering, but not to exclude EMIC waves as a potential 
wave responsible for the proton scattering. The relative contribution 



to the proton aurora emission from MS or EMIC waves deserves a 
future study. 

Finally, in a departure from the previous works^^, we focus on the 
pitch- angle scattering by the MS waves instead of the instability of 
MS waves in this study. The enhanced MS waves appear to be excited 
by those injected anisotropic protons with a typical ring distribution 
at energies of —10 keV or above prior to this event. Unfortunately, 
there is no direct data on MS waves or energetic proton distributions 
before this event because the Cluster satellite doesn't stay in the 
radiation belt. We leave the instability of MS waves to a future study. 

Methods 

The ray tracing of waves is performed by using the following standard Hamiltonian 
equations^*': 

dR _ dD /dD 
'dt~~~dk/d^ 

dk_dD /dD 

where R, co, k, and t represent the position vector of a point on the ray path, the wave 
frequency, the wave vector, and the group time, respectively. The wave dispersion 
relation D obeying D(R, co, k) = 0 at every point along the ray path, has well been 
documented in the previous work. The spatial variation in D can be written: 

dD _ dD dBo dD dN^ dD dk 
'dR~dBo^^dNc^^^dR 

where Bq is a ambient magnetic field and Nc is the background plasma density. 

Here, two Cartesian coordinate systems are adopted for the ray-tracing calcula- 
tion^^'^\ The first is Earth centered Cartesian coordinate system (OXYZ), in which Z 
axis points north along the geomagnetic axis; and the X and Y axes stay in the 
geomagnetic axis equatorial plane. The second is a local Cartesian system, in which 
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Table 1 


Parameters for Ray Tracing 








f(Hz) 


5.0 
5.5 
6.0 


11, 12, 13, 14,23,27,31,35 

12, 14, 16, 18,30, 35,40, 45 
16, 18,20, 22,40, 45,50, 55 


41 .8, 45.5, 49.3, 53.0, 87.2, 1 02.4, 1 1 7.5, 1 32.7 
34.2, 39.8, 45.5, 5 1 .2, 85.4, 99.7, 1 1 3.9, 1 28.2 
35. 1 , 39.5, 43.9, 48.0, 87.8, 98.7, 1 09.7, 1 20.7 



the z axis points along the direction of the ambient magnetic field, the x axis is 
orthogonal to the z axis and stays in the meridian plane pointing away from the Earth 
at the equator, and the y axis completes the right-handed set. The wave vector k makes 
an angle 6 with the z axis and the projection of k onto the xy plane makes an angle rj 
with the X axis, viz., k = k cos 6z-\-k sin 6 cos rjx-\-k sin 6 sin rjy. rj = 0°, 90°, 180° 
and 270° correspond to the perpendicular component k^ pointing away from the 
Earth, toward later MET (eastward), toward the Earth, and toward earlier MET 
(westward), respectively. 

For ray- tracing, we adopt a dipole magnetic field model and the global core plasma 
density modeP^. Considering that MS wave propagate very obliquely, we choose the 
initial wave normal angle 6 = 88°, and the initial azimuthal angle rj = 150° (toward 
earlier MET) and 30° (toward later MET). The other ray- tracing parameters for 
different wave frequencies and source locations are shown in Table 1 . 

To calculate the diffusion rates, we assume that the MS wave spectral density Bj 
follows a typical Gaussian frequency distribution with a center/^, a half width 5f, a 
band between /i and/2^°. 



2B] 



exp 



{f-fn,f 



(4) 



here Bt is the wave amplitude in units of Tesla and erf is the error function. 

We choose the wave normal angle distribution to also satisfy a standard Gaussian 
form: 



|exp[-(X-X^)/Xj^ 



for Xi<X<X2, 
otherwise, 



(5) 



where X = tan 0 (0i < 0 < 02> ^1,2 = tan 0i,2)> with a half- width X^, and a peakX^. 
Based on the observation, we choose X^ = tan 89°, X^o = tan 86°, Xi = X^ — X^j, X2 
= X^ + Xoj; and the maximum latitude for the presence of MS waves 1^ = 5°. We 
assume that the wave spectral intensity remain constant along the dipolar geomag- 
netic field line. 

As shown in Figure 3, there is a cross region between two bands with the cross 
frequency /c;, = 70 Hz. We consider contribution from harmonic resonances up to 
w = ±20 for the lower band and = ±30 for the upper band. To avoid repeating 
calculation from the cross region of two bands, the lower band stops at/^^ = 70 Hz; 
and the upper band starts dXfcr = 70 Hz. We then compute MS-driven bounce- 
averaged diffusion coefficients at the location L = 4.6. 

The evolution of the proton phase space density/ is calculated by solving the 
bounce-averaged pitch angle and momentum diffusion equation 



dt GpdoLe 

1 d 



+ 



Gdp 



1 df 



(Dap) 



G[ {Dp,) 



(6) 



herep is the proton momentum, G = p^T{oLe) sin cos with being the equatorial 
pitch angle, the normalized bounce time T{ae) ~ 1.30 — 0.56 sin a^; (D^a), (Dpp), and 
(Dap) = {Dp^} denote bounce-averaged diffusion coefficients in pitch angle, 
momentum and cross pitch angle- momentum. The explicit expressions of those 
bounce-averaged diffusion coefficients for MS waves can be found in the previous 
work^'*'^^ 

The initial condition is modeled by a kappa-type distribution function of energetic 
protons^^: 



fo{p,sm a,; 



«pr(K; + /+l) 



■ 7^3/203 K;('+V2)r(/+l)r(K;- 1/2) \^ Q 



p sin a 



1 + 



-{K + l+l) 



(7) 



here Hp, I, k, and T denote the number density of energetic protons, the loss-cone 
index, the spectral index, and the gamma function. 6^ represents the effective thermal 
parameter normalized by the proton rest mass energy nipC^ (— 938 MeV). 

Solution of the diffusion equation (6) requires to choose the appropriate and 
realistic initial and boundary conditions in order for a realistic simulation of this 
event. Boundary conditions in pitch angle are taken:/ = 0 at the loss-cone = ocl (sin 
ccl = L~^'^{4 — 3/L)~^'^) to simulate a rapid precipitation of protons inside the loss 
cone (Fig. Id), and df/dcCe = 0 at = 90°. For the energy boundary conditions, /is 
assumed to remain constant at the lower boundary 0.5 keV and the upper boundary 
20 keV, respectively. 

Based on the observational data, we choose the following values of parameters: 
9l = 0.5 X 10"^( ~5 keV), / = 0.01, ?c = 3 and = 2.9 cm''. We solve the equation 
using the recently developed hybrid difference method^^, which is efficient, stable and 



easily parallel programmed. The numerical grid is set to be 91 X 101 and uniform in 
pitch angle and natural logarithm of momentum. 
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